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Abstract 

Herein, we present a canonical form for a natural and necessary generalization of the Lambert W 
function, natural in that it requires minimal mathematical definitions for this generalization, and nec- 
essary in that it provides a means of expressing solutions to a number of physical problems of funda- 
mental nature. In particular, this generalization expresses the exact solutions for general-relativistic 
self-gravitating 2-body and 3-body systems in one spatial and one time dimension. It also expresses 
the solution to a previously unknown mathematical link between the linear gravity problem and the 
quantum mechanical Schrodinger wave equation. 

AMS Numbers: 33E30, 01-01, 01-02 

Also related to: 70B05, 81Q05, 83C47, 1 1A99 



1 Introduction 

The Lambert W function satisfying W{t)e w ^> = t was first introduced by Johann Heinrich Lambert 
(1728-1777), a contemporary of Euler. Though it is more than two hundred years old, its importance 
and universal application was only realized in the last decade of the 20th century. With the combined 
efforts of Gonnet and others, the W function has become a tool. Ironically, it had been injected into 
the Maple Computer Algebra system by Gaston Gonnet over many objections because it did not appear 
as a "standard" special function known in literature (e.g. see [1-3]) though it was useful as a means of 
expressing solutions to transcendental algebraic equations. A presentation of the work of Scott et al. [4] 
in 1992 showed that it expressed an exact solution to a fundamental problem in quantum mechanics. 

This encouraged Corless et al. to make a literature search of the W function to find that it had been 
"invented" and "re-invented" at various moments in history. Its applications were numerous [5]. For 



example, the W function has appeared in electrostatics, statistical mechanics (e.g. [6]), general relativity, 
inflationary cosmology (e.g. [7]), radiative transfer, Wien's Displacement Law of blackbody radiation 
(e.g. [8]), quantum chromodynamics, combinatorial number theory, fuel consumption and population 
growth (see e.g. [9] and references herein) etc. Within the past decade, the Lambert W function has 
been embedded in other computer algebra systems, multiplying the applications of this function and also 
increasing its awareness and thus vindicating its appearance within a Mathematical software system. 

More recently, the Lambert W function has also appeared in the "lineal" gravity two-body problem [10] 
as a solution to the Einstein Field equations in (1 + 1) dimensions. The Lambert W function appears as 
a solution for the case when the two-bodies have exactly the same mass. However, the case of unequal 
masses required a generalization of Lambert's function [10, eq.(81)]. Subsequent discussion brought the 
realization that this generalization for unequal masses had a one-to-one relationship with the problem of 
unequal charges for the quantum mechanical problem of Scott et al. 

This realization fueled the impetus for investigation into a proper generalization, the focus of this article. 
Moreover, it became clear that something vital about the W function had been missed. The information 
and awareness of the literature on Lambert's function is still fragmentary. For example, long before 
D.E.G. Hare [5, (a)] extended the definition of the W function into the complex plane, such an analysis 
had already been done indirectly by Byers Brown [11, 12] not only for the standard W function but also 
for its generalization discussed herein. We see that no matter how exhaustive a literature search is made, 
it cannot address or cover all the aspects of a function over two hundred years old! 

The goal of the present work is to examine this generalization while clarifying a number of issues with 
regards to the unfortunate "fragmentation" of information and awareness of Lambert's function. Of course, 
one can define generalizations in myriad ways. Thus, we seek a generalization that is "natural" i.e. 

1. It is economical in that minimizes the need for new mathematical definitions. 

2. It has applications in nature. Better still, it is ubiquitous to nature, not unlike the standard Lambert 
W function itself. 

3. It expresses solutions to a broad range of mathematical problems. 

4. Its capacity for further generalization and its reduction and correspondence to the standard Lambert 
W function is transparent. 

Such a function satisfying these criteria is clearly a fundamental mathematical structure worthy of consid- 
eration in the Mathematical/Physical literature. 

In this work, we present a canonical form for a generalization of the W function which satisfies this 
criteria. This is done as follows. First, we then re-examine the quantum mechanical problem of Scott et 
al. whose solution also expresses the solution to the linear gravity problem of Mann and Ohta [10]. We 
solve the case of unequal charges (unequal masses for the gravitational problem) intuitively. The impetus 
is partially derived from the notion of P.A.M. Dirac that a sound mathematical structure has a potential 
basis in reality (and the converse might just be true!). 

Next we find that the generalization fits into a tetration framework (or iterative exponentiation), and re- 
quires only a nesting of the definitions of the Lambert W function, thus satisfying the first requirement. 
Next, we seek solutions to the gravitational three-body problem in (1 + 1) dimensions and find that the 
first generalization can be naturally extended further. Finally, it is found that the end result expresses 
solutions to a huge class of delayed differential equations. It is also helpful in expressing the solution to 
the (three-dimensional) hydrogen molecular ion i.e. the quantum-mechanical three-body problem for the 
case of clamped nuclei. Concluding comments are made at the end. 



2 One-Dimensional Quantum Problem 



The one-dimensional version of the hydrogen molecular ion [11, 12, 14] is given by the double Dirac 
delta function model: ^ 

-\^- q ^ + X5{ ^ x - RM = E(\W (i) 

where Za = q and Zb = A q. The ansatz for the solution has been known since the work of Frost [15]: 

$ = Ae~ d \ x \ + Be~ dlx - Rl (2) 

where < R < oo. All quantities are real. Matching of ij) at the peaks of the Dirac delta functions at 
x = 0, R yields: 

q — d qe~ dR 
q \ e -dR q X-d 



(3) 



and the energies are given by E± = —d\/2 where d± is governed by the secular determinant of eq. © 
when it is set to equal zero ( [4, eq.(17)] for q = 1): 

d±(A) = \q{\ + l)±\{q\l + \?-A\q 2 [l-e-^ R ]} 1/2 (4) 

When A = 1, the pseudo-quadratic in ((4) reduces to: 

d± = q[l ± e~ d±R ] (5) 

Although, the above has been known for more about half a century, it was not until the work Scott et al. [4] 
that the solution for d± was exactly found to be: 

d± = q + W(0,±qRe- qR )/R (6) 

where ± represent respectively the symmetric or gerade solution and the anti-symmetric or ungerade 
solution. The first argument of the W function, being zero, reminds us that Lambert's function has an 
infinite number of branches and that we are selecting the principal branch. 

The anti-symmetric case, cL is interesting because it appears to go to zero as R — > 1 for 
q = 1 , in other words the energy goes to zero and the corresponding eigenstate appears to go into 
the continuum. However, for R < 0, W{— 1, —Re~ R ) which has an order 2 branch at R = 1 yields a 
real number. Given the analysis in the complex plane of the energy eigenstates [11, 12], we can see there 
was already awareness of more than one branch for the solution as far back as the 1970s by mathematical 
physicists well versed in the mathematical framework of linear molecules. 

So far, the Lambert W function could only express the solution for the case of equal charges. We now 
examine the general case of unequal charges. The pseudo-quadratic of eq. (@) seems complicated until we 
rewrite it in a very simple form (for q = 1): 

_ 2rR (1 — x) (X — x) 

e 2xR = ~y where x = d± . (7) 

A 

The above encapsulates both "gerade " and "ungerade " solutions 2 . It also represents a canonical form for 
a whole class of transcendental algebraic equations. The right side of © is a quadratic polynomial in x 
only while the left-hand side is a function of x and R. It must be emphasized that R is a constant in the 
range [0, oo) consequently allowing infinite choices for R. The left-hand side of CZ) is a whole parameter 
family of curves in x while the right-hand side represents only one curve in x. Therefore, eq. (0 is a 

'For simplicity, we set the charge q — 1 for the rest of this work. 

"Symmetry is lost when A / 1 and the terms "gerade" and "ungerade" no longer have the same meaning. 



canonical form for an implicit equation for x. Note that when A = 1, we have a double root for this 
polynomial and both sides of (f7) factors into two possible cases where the solutions are given by eq. © 
for q = 1. The problem in linear gravity [10] namely eq. (82) of ref. [10]: 

y 2 = a 2 + {x 2 — a 2 ) exp (2x) exp (— 2y) 

relates exactly to eq. (£7) by the following transformation: 

2x x — y 
A = 1 where R = — (x + a) and d = (8) 

x + a x + a 

when q = 1. The case A / 1 represents the generalization we seek. Thus, we seek a solution to: 

e~ 2xR = a b (x - n)(x - r 2 ) (9) 

where in relation to the above problems {r±, r 2 } = {1, A} are the real roots of a quadratic polynomial 
and where {a , b } = {1, 1/A}. However, we treat these parameters generally while making all necessary 
assumptions to ensure a real solution. Since a quadratic is merely a product of first order polynomials, this 
guides us intuitively to consider the following. We assume there exists a value y such that: 

e'* 3 * = a (x-n) (10) 

e -Rx(2-y) = 6o ( x _ r2 ) (11) 

Multiplication of the left sides and right sides of eqs.([T01 and (IT)) yields ([9) the equation we desire to 
solve. However, individually eqs.CTOl and ([ID can be solved using the standard Lambert W function: 

W(yRe-^y R /a ) 

a (x 1 -r 1 ) = a a (12) 

yR 

, , W{{2-y)Re-^ 2 -y) R /b ) 
b (x 2 -r 2 ) = bo ^— ^ (13) 

Letting y = 1 + e, we seek y such that x\ = %i. Thus, the "separation" parameter is governed by: 

. W((l-e)Re- r ^ l -^ R /b ) Will + e) Re~ r ^ 1+ ^ R / a ) 

(ri " r2) = o^7)R (TT7)R (14) 

Substituting eqs.(12) and (B) into © 

_ 2xR _ W ((1 + e)Re- r ^ R /a ) W ((1 - e)Re- r ^ R /b ) 
6 " a ° b ° (l + e)(l-e)i?2 {n) 

and taking logarithms on both sides of ([15) allows us to isolate an expression for x subject to the constraint 
that y is such that x\ = x 2 . Looking at "canonical" forms of eqs.(T2) and CT3) in comparison with ( [T5) 
makes us infer the generalized Lambert W function as 

Q 2 = n 2 (a ,b , ri ,r 2 ,R) = W{ Zl )W{z 2 ) (16) 

where 

Zl = (l + e)Re- r ^ 1+ ^ R /a 
z 2 = (l-e)Re- r ^~^ R /b 

and where e = e(a , b ,r\,r 2 ,R). The above is a product of standard Lambert W functions in the same 
fashion a quadratic polynomial is the product of first order polynomials. When r% = r 2 , it is clear that 
e = (or when an analytical expression for e is possible) and we recover the solution in terms of the 
standard Lambert W function. The subscript 2 on Q reminds us that the right side is a second order 



polynomial. In general, y or for that matter e will be referred to as a separation parameter, which allows 
the generalized function shown here to be decoupled as a product of standard Lambert W functions. Note 
that when a = b = 1, we recover the symmetric (gerade) solution and with a a = b a = —1, the 
anti-symmetric (ungerade) solution of eq. @. 

In view of eq. (141) , one realizes that the separation parameter e is itself governed by a transcendental 
equation which looks even more complicated than the original transcendental algebraic equation of eq. 0. 
The critic then asks: how are we further ahead? Part of the answer lies in considering an important aspect 
about generalizations i.e. whether or not they have a capacity to collapse into special cases other than 
the original function from which the generalization was inferred or collapse into previously unknown 
solutions. Let a Q , b a ,ri, T2 be the values as quoted below ([9]) and let us make a series expansion of x\ — X2 
in the parameter A: 

1 + e W(R{l+e)e- R ^+*)) R(l+e) 

xi-x 2 ~ 7 stt-ts 2A - 2R(-1 + e) A 2 

28 

-4# 2 (-l + e) 2 A 3 - -^i2 3 (-l + e) 3 A 4 -24# 4 (-l + e) 4 A 5 + 0(A 6 ) (17) 
By inspection, we can see that e = 1 eliminates all terms of order A greater than 1 and consequently: 

1 +e W(2Re- 2R ) e 2R 
e W(2Re~ 2R ) e 2R 

Solving for A such that the difference x\ — X2 is zero, yields: 



xi-x 2 = „„„ p „_ 2fl x no 2A 



1 W(2Re-'«) 

2 4R 



In this case, the solution for x is found to be: 



1 , fW(2Re- 2R ) 

In — 

2R V 2R 



In 1 v , n '-) (19) 



Thus for A satisfying (TT8T) and e — > 1, we have a previously unknown solution to eq. (7) in terms of the 
standard Lambert W function. Granted e — ► 1 represents a limiting extreme: this solution has been vin- 
dicated by numerical and analytical demonstrations using computer algebra. Thus, we have a previously 
unknown particular solution for the case of unequal charges (or unequal masses in the linear gravity prob- 
lem) but for a peculiar value of A dependent upon R which admittedly is not physically useful since the 
physical charges Za and Zb are constants independent of R. This bears a striking resemblance of the 
results of Demkov who found analytical solutions for the three-dimensional hydrogen molecular ion 
but for a particular choice of charges (again not physically useful for the same reasons) which reduced to 
Whittaker functions i.e. the type of solutions found for the hydrogen atom [16]. 

At any rate, this demonstration shows that the generalization we inferred is not impotent but can also 
collapse into simpler special functions for special cases. There are other such cases. Since eq.([14l has 
the form r\ — r2 = ffa) — <?(n), one particular obvious solution is simply r\ = f{r2) and T2 = g(r\). 
Decoupling this particular solution with respect to n and T2 yields equations governing the latter: 

W{R{1 - e) exp( * L )/b ) 

n = W^e) (20) 



T2 



WW. + e) exp( (1 -^( (1 -t; 2ii(1 ' £)/fco) )/a ) 



R(l + e) 

At this stage, we can change our context: rather than look for solutions of eq. ® for arbitrary choices of 
the parameters R, a a , b a , r\ and T2 and be burdened with the analytical determination of e: we instead use 



e as a common parameter to r% and r 2 allowing us to find which values of r\ and r 2 satisfy eq. (9) for a 

given choice of the remaining parameters a Q b Q and R. Moreover, we also find the following previously 
unknown exact solutions to ®: 

-2fl(l + ft r 2 ) 

1 e b o l + b Q r 2 

r i = 7-, a Q = , x = , e->l. (21) 

o r 2 b a 

where r 2 and b a are arbitrary real numbers and also 

T2 = h a ° = -2r lR lHr 1 b o y I= ^ k(riW ' (22) 

where n and b a are arbitrary real numbers. Eq. ( [21) is in closed form (elementary functions) and ([221 
is in terms of the Lambert W function which bears some resemblance to eq. (T9T ) and could be called 
another "Demkov"-type solution. This is one of the utilities of the "separation" parameter e: eqs. (14]) and 
especially (20) look more complicated than the original equation (9), they are however useful in finding its 
particular solutions 3 . The parameter e also illustrates the reduction from the proposed generalized function 
to the simpler (standard) Lambert W function in a transparent manner. In the next section, we properly 
define the function Q, 2 which we derived intuitively. 



3 Iterated Exponentiation 

Infinitely iterated exponentiation or tetration is defined as the limit, 

lim a n = a° a ' = Sj(a) (23) 

n^oo 

This function can be written compactly as: 

Sj(a) =e~ w{ - lna) aeC (24) 
The product on n W functions is given by: 

W( Zl )W(z 2 ) . . . W(z n ) = (z lZ2 . . . Zn ) e -[W(z 1 )+W(z 2 )+...+w(z n )] 



following from the defining relation z = W(z)e w ( z > . The Lambert W function is governed by the addition 
law [17]: 

W(a) + W(b) = W(ab [l/W(a) + l/W{b)\) (25) 
for K(a) or > 0. By induction, the product is then: 

W( Zl )W(z 2 )...W(z n ) = ( Zl z 2 ...z n )e- W ^ 

= (z lZ2 ...z n )W(f n )/f n (26) 

where stands for an expression involving W functions of z\, z 2 ■ ■ ■ z n that follows from the addition 
law. E.g. for n = 2 and n = 3, 

h = z lZ2 (l/W( Zl ) + l/W(z 2 )) 

h = f2Z3(l/W(f 2 ) + l/W(z 3 )) 



3 We anticipate the existence of other cases when e can be solved in closed form but of course, in general, this is not the case. 
Indeed, e might not even exist which is why the generalization is needed. 



Thus, 



W( Zl )W(z 2 ) = z 1 z 2 W(f 2 )/(z 1 z 2 (l/W(z 1 ) + l/W{z 2 ))) 

= W(f 2 )/(l/W( Zl ) + 1/W(z 2 )). (27) 

In the context of the tetration function Sj(a) of ([24]), the product of Lambert W functions is really just: 

n n (z lt Z2, ...,z n ) = W( Zl )W(z 2 ) . . . W(z n ) = {z x z 2 . . . z n )fi(e- f ") (28) 

a multiple of the continued tetration function. To reiterate, for eqs.(fl2) and (13) when x = x\ = x 2 , the 
solution to eq. (§) is given by a product of W functions, as one substitutes Xj into (x — rj) for j = 1,2. 
Hence, the solution is in terms of one W function, albeit nested. In this framework, thanks to the addition 
law for the W function, no generalization in the sense of additional mathematical definitions is needed: it 
is simply one single W function evaluated at a point involving other W functions. 

Furthermore, if we use eq. ( [28] ) for n = 2 and combine with the earlier result in (I6[) , we obtain a simple 
yet general relation between the separation parameter e and x: 

(r 2 - n)e = + (n + r 2 ) - 2 x (29) 

which is consistent with eq. (15) in view of the addition theorem for W functions. We can already identify 
regimes. For example, for real positive roots T\ and r 2 the term W(f 2 ) is bounded in many cases and 
consequently 4 liniR^oo e = ±1. 

Although we initially assumed 3ft z n is not a negative number, injecting of numbers into eq. ([281 show that 
this identity formula is very robust after all if one considers different W functions on different branches. 
For example, ([28]) holds if x\ = 1 + 2i, x 2 = —1 for W(x 2 ) = W(— 1, x 2 ) and is complex-valued. 



4 Rational Polynomials 



In the previous sections, we considered the right side of the transcendental equations to be polynomials, 
but one can also consider a rational polynomial. Let us consider the equation where the right side is a ratio 
of first order polynomials: 

e -2R* _ 

In parallel to what we did before, we can consider 

e -Rx{2-y) 

Eq.([32) can also be solved in terms of the standard Lambert W function i.e. once the inverse is taken on 
both sides of ([321 i.e. (2 — y) — > (y — 2) in eq. ([IT). The generalization is thus: 

= W ((1 + e)Re- ri ^ R /a ) W (-(1 - e)i?e +r2 ( 1 - e ) ii /6 ) (33) 

where the subscripts 1,-1 denote respectively the polynomial degrees in ( [30T ). This also fits into the 
tetration framework of section [3J the only change being done to the argument x 2 and the right side of (30[ ) 
can be generalized to a rational polynomial of higher order as both numerator and denominator can be 
"separated" in a manner demonstrated in eqs.(3T]) and ( [321 to yield generalized Q functions that can be 
further nested together using the addition theorem. The next section presents applications of this type of 
problem. 

4 Of course, we realize that it is much easier to solve the transcendental equations numerically; the goal of this exercise is to 
define and justify our generalization of the W function. 



a (x - ri) 
b (x - r 2 ) 



a Q (x - n) 
1 

b (x - r 2 ) 



(30) 



(31) 
(32) 



5 Three-Body Linear Gravitational Motion 



The solution of the three-body one (1 + 1) dimensions via dilation theory requires solving for V which is 
governed by the following equation [18, eq.(32)]: 

(V - mi) (V - m a ) (V - m 3 ) (34) 
= mi m 2 (V — s% s 2 777,3) exp(K V R \ sm(q)\) 

( TT 

+ mi m^(V + sq s 2 m 2 ) exp(iv~ V R | sin ( q + — 



7T 

+ 1712 ^3 (V — sq si mi) exp(K V R | sin ( q — — 



(35) 



where 



si = sign(sin(<7 + |)) 
s 2 = - sign (sin(g - |)) 
sq = sign (sin(g)) 

This problem has no closed form solution. The hard part is in trying to extricate the exponential terms 
from the rest of the expression. At q = or ±7r/3, the trigonometric quantities s, for i = 1,2 and sq 
simplify to specific values. E.g. for q = and V 7^ 0, it is found that: 

exp(-2R t V) = — ?- r (V - m 3 )(F - (mj + m 2 )) (36) 

m 3 (mi + m 2 ) 

where i? 4 = KRy/l/4, As we can see V is governed by exactly the same type of transcendental equation 
as eq. ® whose solution, as we have seen, can be expressed in terms of our f2 2 function. In the case when 
7773 = mi + m 2 , we have a double root and thus a solution for 06] ) can be expressed exactly in terms of 
the standard W function and has two solutions: 

V = 1TI3 — —W (± Rt m^e Rtm:i ) when q = 0, m 3 = mi + m 2 (37) 
Rt 

Very similar solutions also exists for q = tt/3 with mi = m 2 + m 3 and for q = —tt/3 and q = 2tt/3 with 
m 2 = mi + m 3 . This demonstrates that at every tt/3, solutions in terms of our generalized Q functions 
exist. At q = tt/6 and equal masses i.e. mi = m 2 = m 3 we arrive at two possible equations depending 
on the outcome of factorization. One equation is: 

(V - m) 

exp(cV) = >- , (38) 

m 



and the other equation is: 



(V — m) 2 

eMcV) = - { > (39) 

m{V + m) 



where c = KR/2. Eq.([38) can be solved in terms of the standard W function i.e. 

W(mce m 



V = m 1 



mc 



which is in the same form as the range x range and drift Zdrfft equations of section 2 and the anti-symmetric 
solution d- in section 3 with the exact correspondence given by rj = —mc and b = f Xj o- However the 
right side of 091 involves a rational polynomial of the form P 2 (V)/Qi(V~) and requires fi 2; _i to express 
the solution. If we consider the region < q < tt/3, we obtain the form: 

exp(-c v/(3) cos(g)y) = where M,N^oo (40) 

Qm{V) 



This rational polynomial can be generated from: 



fjy- mg) (V2-V* K + m 2 ) - mi * m 2 {f 2 - 1)) ^ P N (V) 
777,3 * (mi * ni2 * (f 2 — 1) + mi * / 2 + m-2 * V) Qm(V) 

where / = exp(c sin(q)V) generates the rational polynomials. In this regime sin(g) q is small and 
cos(q) ~ 1 and thus this treatment is possible. A similar story applies to other regimes [vr/3, 2tt/3], 
[27r/3, 7r] etc . . . Thus, we can see that the general solution falls in the following general form. 



6 General Form 



Thus, the fully generalized form concerns expressing solutions to this general class of transcendental 
algebraic equation: 

Qm[x) 

where k > is a constant and Pjy(x) and Qm{x) are polynomials in x of respectively orders N and M. 
The general solution can be expressed by which is formally a product of N+M (standard) Lambert 

W functions with N + M — 1 "separation" constants. 

The standard W function applies for cases when N = 1 and M = and expresses solutions for the case 
of equal charges for eq. CD or the case of equal masses for the lineal two-body (1 + 1) gravity problem. 
Correspondingly, the case N = 2 and M = expresses solutions for the case of unequal charges (or 
unequal masses for the lineal two-body gravity problem) and some particular cases of the lineal three- 
body gravity problem. In the limit as M, N — > oo, this equation can be used to express solutions of the 
three-body lineal gravity problem as we have just shown. 

Moreover, the case N = 2 and M = (and more generally N = 2 and M = 1 ) express the solutions of 
a significant class of delayed differential equations [19, eq.(3)]. These arise in a variety of mechanical or 
neuro-mechanical (oscillatory) systems in which non-linear feedback plays an important role. These have 
applications in e.g. models for physiological systems (medicine) [20]. 

Recently Adilet Imambekov and Eugene Demler [21] considered a bose-fermi mixture in one dimension. 
We note that their Eq.(2): 

N d 2 

H = -Y.q^2 +2cJ2^ i -x j ) c>0 

i % i<j 

is in fact the Schrodinger wave equation for a linear molecule. Their analytical solution as written in 
their eq. (20) is indeed a special case of the general form expressed (14TT) . Thus, we can also see that this 
function has fundamental applications in Physics and can play a fundamental role in Mathematics. 



7 Conclusions 



Thus we have identified a generalization of the Lambert W function or Omega function as solutions to a 
large class of transcendental equations as written in eq. ( |4TT) . This generalization, denoted Q n , satisfies 
the criteria mentioned earlier in our introduction by using the analytic framework of tetration. 

We have also shown that the two-body problem in lineal gravity and double well linear quantum mechanics 
have the same generalization of W, namely &2 as solutions. The reason why is because the linear gravity 



theory via dilaton theory produces a partial differential equation, namely eq. (30) of ref. [10] which can be 
treated formally as the Schrodinger wave equation as written in (JT). This is to be elaborated elsewhere [22]. 

Furthermore, this work on the Lambert W function has helped in finding analytic solutions to the quantum 
mechanical 3-body problem known as the hydrogen molecular ion in the case of clamped nuclei of equal 
charges [23]. All this vindicates our proposed generalization of the W function as being of fundamental 
and physical importance. 

Most of the special functions in the known literature (e.g. [1]) are special cases of the hypergeometric 
functions and/or the Meijer G-function [3]. The Lambert W function apparently bears no relationship to 
these functions and belongs to a class of its own. The generalization we have presented is a first step in 
identifying that class. 

Although we have inferred a canonical form for a generalization as expressed by ( [41]) and given mathemat- 
ical and physical justifications for it, we have yet to clearly identify a domain and range of applicability or 
conditions of analytic continuation. Neither have we formulated Taylor series nor asymptotic series useful 
for computation. Naturally, it is expected that our "separation" parameter e in (J5) likely has a restricted 
domain of applicability (its primary use being to infer the general result starting from the standard W 
function). 

Nonetheless, given that we have fast computational means to transcendental algebraic equations, the task 
of obtaining the floating-point attributes in many cases is trivial today with modern algorithms. Ironically, 
it has been our capacity to readily solve these transcendental equations numerically that has made many 
take an analytical solution for granted. This may have been a mistake as we can see that such analytical 
solutions are ubiquitous to certain fundamental problems in Physics and Mathematics. In a true sense, this 
work is only a beginning. 
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